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1 Introduction 

In this paper we present an error analysis for some Newton-Cotes quadrature 
formulae. We consider Simpson's rule, 3/8 Simpson rule and Boole's rule. A 
similar error analysis for Simpson's rule have been investigated more recently 
(D> El EH j HI; El) w ith the view of obtaining bounds on the quadrature rule 
in terms of a variety of norms involving, at most, the first derivative. It is 
well known that if the mapping / is neither four times differentiable nor is the 
fourth derivative bounded, then we cannot apply the classical Simpson's 
quadrature formula, which, actually is one of the most used quadrature formulas 
in practical applications. Thus, the above mentioned analysis is important as 
well as the analysis presented here. 

The current work brings results for the above mentioned Newton-Cotes 
quadrature rules giving explicit error bounds and using results from the modern 
theory of inequalities. The used inequalities are known in the literature as in- 
equalities of Ostrowski-Gruss type. The error bounds are expressed in terms of 
second derivatives. As we have already mentioned for Simpson's rule, the gen- 
eral approach used in the past involves the assumption of bounded derivatives 
of degree higher than two. We also mention that the obtained results can be 
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derived using the Peano kernel theorem. In any case, these bounds are gener- 
ally, but not always, better that the usual Peano error bounds (see Remarks [5] 
□ and©. 

Here we do not consider composite quadrature rules since they can be formed 
in the usual way. However, the analysis presented here allows the determination 
of the partition required that would assure the accuracy the result would be 
within a prescribed error tolerance. 

In Section 2 we established some auxiliary results. We use these results 
in further sections. In Section 3 we consider Simpson's rule. In Section 4 we 
consider 3/8 Simpson rule and in Section 5 we consider Boole's rule. 

2 Preliminary results 

Lemma 1 Let I C R be an open interval and a, b E /, a < b. Let f : I — ► R 

be a twice differ entiable function and let x E [a, b] be a fixed element. Then we 
have 

b 

f(x)(b - a) - (x - i±*) [/(&) - /(a)] - / f{t)dt (1) 



b b 

1 



b — a 



p(x, t)p(t, s)f"(s)dsdt, 



a a 



where 



Proof. Integrating by parts, we have 



b b 

p(x,t)f'(t)dt = f(x)(b-a)- f f(t)dt, 



l. e. 

b b 

f(x) = ^i- J f(t)dt + ^i- J P (x,t)f(t)dt. 

a a 

If we substitute / — > /' in the above relation, then we get 

b 

fit) = {f(b) /(a)] + ^ I p(t,s)f"(s)ds. 



2 



Thus, 





J p(x,t)f'(t)dt 



p(x,t) 



' [f(b) - f(a)} + ^- I ,,U.s)f" ( s),ls 



b — a 



dt. 



We also have 



I SI /I \/ a + ^\ 

p(x,t)dt = (b-a)(x — ). 



From the above relations it follows 

b b 



1 



b- 



p{x,t)p(t,s)f"(s)dsdt 



o 

f(x)(b -a)-{x- i±^) [/(&) - f(a)} J f(t)dt. 



Lemma 2 Letp(x,t) be defined by 0). Then we have 
b 

a 

f (6-a)(a;-^)( S -a)-^( S -a) 2 , se[a,x] 
{ {b-a)(x-^)(s-b)-!^i{s~b) 2 , se(x,b] 

where x £ [a, 6] . 

Proof. For s e [a, x] we have 



(3) 



q(x,s) 



(t -a)(s- b)dt + (t- a)(s - a)dt + (t - b)(s - a)dt 



. ,.(s-a) 2 , (x - a) 2 - (s - a) 2 , (x - b) 2 

(s - b y - + (s - a)± 1- - ( s _ a ) y - 

a+b b-a 2 
(6 - a)(x — )(s - a) — (s - a) . 
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For s G (x, b] we have 



g(x,s) 



(f - a)(s - b)dt + (t- b)(s - b)dt + (t-b)(s- a)dt 



Ax -a) 2 , L .(x - b) 2 + (a - 6) 2 , . (a - 6) 2 

- (6-a)(x- -^-)( fl _6) - __( s _6) 2 . 
From the above relations we see that (j3J holds. ■ 

Corollary 3 Let the assumptions of Lemma\^ be satisfied. Then we have 



q(a,s) = 
q(b,s) 



II ■ h 
,3a + b 



a + 36 . 



,2a + 6 , 

q{ — o — > s ) 



,a + 2b 

q{ — o — > s ) 



-\{b-a)*(8-b)-^{a-b)*, 

1 - {h -af{s-a)- b -^( S -a)\ 
- b -^{s-a) 2 , S e[a,f 



b—a 



(s-b) 2 , se{^,b] ' 



3a+&1 



-\{b - a) 2 (s - a) - ^(s - a) 2 , s G [a, 
^(&-a) 2 ( S -&)-^( S -&) 2 , ae(^>J ' 

\{b - a) 2 (s - a) - - a) 2 , s G [a, ^1 

±(6-a) 2 (s-&)-^(s-6) 2 , se(2f,t] ' 

-i(6- a ) 2 (.s-a)-^( S -a) 2 , sek^J 
-i(6- a ) 2 ( S -6)-^( S -6) 2 , *G ' 

|(6 - a) 2 (s - a) - *=2(a - a) 2 , a G [a, ^] 



3 Simpson's rule 

Theorem 4 Let the assumptions of Lemma^hold and let 7, T be real numbers 
such that 7 < /"(f) < r, Vi G [a, 6]. TTien u>e /lave 



G 



/(«)+4/(^) + /(6) 



<- -(b-a) 3 . (4) 

" 162 1 ; W 



Proof. For x = a the left-hand side in |Q is equal to 
(6-a)M+M_/ /(t)& 



(5) 
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For x = b the left-hand side in (JIJ is equal to 



(b-a) 



f(b)+f(a) 



f(t)dt. 



For x = the left-hand side in Q is equal to 



If we now multiply UJ by 2 and add (JSJ or © then we get 



b — a 



/(«) + 4/(^) + /(6) 



- 3 J f(t)dt. 



The corresponding right-hand side is 



b b 



R{a,b) = 



1 



b — a 
1 

b — a 



, a I b . , . 
2p(— ,t)+p(o,t) 



a a 
6 



2q(— ^— ,s) +<?(a,s) 



p{t,s)f"(s)dtds 
"(s)ds. 



Let if x (s) = 2g(2±^, s) + g(a, s). Then we have 





y Ki(s)ds = 



such that 



K 1 (s)f"{ a )ds = / Kx(s) 



f"(s) 



r + 7 " 



ds 



and 



\R(a,b)\ < 



1 

max 

b — ase[a,b] 



f"(s) 



r + 7 



\Kx(s)\ds 



since 



max 

s£[a,b] 



/"(*) 



r + 7 



< 



r- 7 



Hence, 



\R(a,b)\ < 



r-7 

2(6 - a) 



2g( — — ,s) + g(a,s) 



ds. 



(10) 



(11) 



We now calculate 





I 



„ .a + 6 

2g(— ^— ,s) + g(a, s) 



fi.s 



(12) 



. 6 — a 
~~2 



{b-af 



b — a 



+ 





/ 



, . 9 (6 — a) 2 , , . b — a . 

-2— (* - &) 2 - ^r^(s - &) - — (s - fo) 2 



.6 — a , 



From the equation 



, 9 (6 — a) 2 , , , b — a . 
-2^— (s - a) - - — ^-(s - b) — (s - b) 2 = 



.6 — a, 
~2~ 



we find the solutions 



From the equation 



b — a 
~~2 



we find the solutions 



Si = a, s 2 = 



2a + 6 



(13) 



{b-af 



b — a 



-2 bf ^-^(s -b)- U —^( S -bf = 



2 x ' 2 
a + 26 



s 3 = 6, s 4 = 



(14) 



6 



From JT2Jl-lO we have 



2q(—^—,s) +q(a,s) 



ds 



b — a 



(b-af 



b — a 



-2—1 ( s _ a f _ A^-^( s - b )- U —^(s- bf 



ds 



+ 



„6 — a, . 9 (b — a) 2 . ,. b — a, , ., 



2 m ■ I 
3 



+ 



2 

6 



6- 



(6 -a) 2 



(6 — a) 2 . , . b — a , > 



.b — a , 
~2~ 



108^ ; 27 v ; 
From the above relation and (f5fl — (| 1 1|> we get 



b — a 



f(a) + if(^) + f(b) 



From i|15|) we easily get Q. ■ 

Remark 5 T/ie usual Peano error bound is 



3 / f(t)dt 



(15) 



b — a 



/(a)+4/(i±^)+/(6) 



/(i)dt 



_ 81 lu 



(16) 



J/ me choose 



7= inf /», r= sup/"(s) 
»e[a,&] se[o,6] 



£/ien < ll/"^ and it is obvious that ^ is better than Mb]) . In fact, these 
two bounds are equal if and only ifT= —7. This case (V = —^) is very rare in 
practice. Specially, ifT is large and V « 7 then ^j) is much better than Mb)) . 
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4 3/8 Simpson rule 

Theorem 6 Under the assumptions of Theorem^we have 



b — a 



/(a) + 3/(^) + 3/(i±*)+/(6) 



f{t)dt 



384 



(17) 



Proof. If we substitute 
left-hand side 



_ 2a+b 



in Q then we get the corresponding 



b -^[f(b)-f(a)] + f(^±^)(b-a)- I /(/),//. 



(18) 



For x = a t? b we have the corresponding left-hand side 



-^[/(6)_/(o)]+/(^±^)(6-o)- //l/.k//. (10) 



If we now multiply (|18|l and l|19|l by 3 and add J5J and (|BJ) then we get the 
left-hand side of the form 



(b-a) 



/(a)+3/(^) + 3/(^)+/(6) 



- 8 / /(t)dt. (20) 



The corresponding right-hand side is 

b 

1 



b-. 



2a + b a + 2b 

q(a, s) + 3q( — - — , s) + 3q{ — - — , s) + q{b, s) 



f"{s)ds. 



Let K 2 (s) = q{a, s) + 3q(^d±, s) + 3q(^,s) + q{b, s). Then we have 

6 



(21) 



J K 2 (s)ds = 



such that 



b — a 



K 2 (s)f"(s)ds 



b — a 



K 2 {s) 



f"(s) 



r + 7 



ds 



8 



and 



1 



\R{a, b)\ < max 

6 — asG[a,b] 



f"(s) 



r + 7 



\K 2 (s)\ds 



< 



r-7 

2(6 - a) 



\K 2 (s)\ds, 



since (|10l) holds. We now calculate 
6 



, , ,2a + 6 . .a + 26 
q(a, s) + 3q( — - — , s) + 3g( — - — , s) + g(6, s) 



ds 



1(6 - a) 3 - 1(6 - a)(* - 6) 2 - Z(6 - a) (a - a) 2 



+ / \2(b-a)([(s-a) 2 + (s-b) 2 ]-(b-a) 3 \ds 



1(6 - a) 3 - 1(6 - a)(* - a) 2 - l{b a){s bf 



2n | /. 

3 

b 



ds. 



From the equation 



1(6 - a) 3 - 1(6 - a)(* - 6) 2 - \{b a)(s - af = 0. 



we find the solutions 



3a + 6 



Si = a, S2 



From the equation 

2(6-a)([(s-a) 2 + (s-6) 2 ] - (6 - a) 3 = 
we find the solution 



S3 



a + 6 



From the equation 
1 



7, 



(6 - ay - -(6 -a)(s- a) 2 - -(6 - a)(s - 6) 2 = 



we find the solutions 



s 4 = b, s 5 = 



2 

a + 36 



From the above relations we get 

6 



„ ,2a + b a + 2b . ,, 

q(a, s) + 3g( — - — , s) + 3g( — - — , s) + g(6, s) 



ds 



\{b - af \(b - a)(s bf 7 -(b -a)(s- af 



ds 



l(b~af-hb-a)(s-bf- 7 -(b-a)(s-a) 



ds 



3a + b 
4 

a + 2b 



+ / [2(6 - a)((s - bf + (s - a) 2 ) - (b - af] ds 



.(b - af - -(b - a)(s - af - -(b - a)(s - bf 



2n ■ /, 

3 

a + 3b 



ds 



i ■ 21. 

3 



■(b - af - -{b - a)(s - af - -(b - a)(s - bf 



ds 



From the above relation and H2()(l and (|22|l we get 



(b -a) 



/(«)+3/(^) + 3/(i±^) + /(6) 



" 48 1 ; 



This completes the proof. ■ 

Remark 7 77ie usual Peano error bound is 



/(-)+3/(H?±*) + 3/(^) + /(6) 



/(t)dt 



II f'll 

<i!^(6-a)». 



-For i/ie reasons given in Remark\^the estimation obtained in Theorem^ which 
is a Peano-like bound, is better than the above Peano bound. 
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5 Boole's rule 

Theorem 8 Under the assumptions of Theorem\^we have 

[7/(o) + 32/(^±&) + 12/(s±&) + 32/(^6) + 7/(6)] 



Jf(t)dt 



<27lff5(r-7)(fc- a ) : 



(26) 



.< = zzft- . We 



Proof. We first write left-hand sides of (JTJ for x — and " — 

have 

6 

f( 3a + b ) + b_a m f{a)] I fm 



(27) 



and 



(28) 



If we now multiply 10 and © by 7, 10 by 12, and (gSJ by 32 and sum the 
obtained results then we get the left-hand side of the form 



(6-0) 



7/(a) + 32/( *!±») + 12/(£±») + 32/(£±») + 7/(6) 



For the corresponding right-hand side we get 
R{a,b) 



-90 / f(t)dt. 
(29) 



(30) 



1 



- / n /3a + 6 , „ ,a + b . , a + 36 . 

7q(a, s) + 32g( — ^ — , s) + 12g(—, fl) + 32<?(^^, s) 



6 — a 
+7#,s)]/"(s)ds. 

Let if 3 (s) = 7g(a, a) + 32q(^, s) + 12 9 (*±*, s) + 32 g (^, s) + 7<z(6, s). 
Then we have 

6 



such that 



b- 



K 3 (s)f"(s)ds = / tf 3 (s) 



/» - 



r + 7 



11 



and 



\R(a,b)\ < 



1 



< 



6 — ase[a,b] 

b 

r-7 

2(6 -a) 



f"(s) 



r + 7 



\K 3 (s)\ds 



\K 3 (s)\ds, 



since (f 101) holds. Hence, 



\R(a,b)\ < 



r- 7 



2(6 -a 

a 

+7q(b,s)\ds 

We now calculate 
b 



« / s o« /3a + 6 ,a + 6 . „„ .a + 36 , 

7g(a, s) + 32g(^^, s) + 12<?(-^-, fl) + 32g(^^, s) 



/ 



,30 + 6 ,a + 6 ,a + 36 . ,, 

7q(a, s) + 32g(^^, s) + 12<?( — , s) + 32g(^^, s) + 7g(6, s) 



ds 



7(6 - a) 



[(« - 6) 2 + (s - a) 2 ] - 38(6 - a)(s - a) 2 + ^(6 - a) E 



a + b 
2 

+ / 


39(6- 


3a + b 
4 




a + 3b 
4 




+ / 


39(6- 


a + t, 
2 




6 

+ / 


7(6 


a + 3b 





ds 



{b-af 



{b-af 



, a + b. 
+ (s - -=-) 



, a + 6, 

+ (* - -s-) 



93 

y(6-a) 3 + 6(6-a)(s-a) 2 



93 

__ (6 _ a) 3 + 6(6 _ a)(s _ 6r 



ds 



[(« - 6) 2 + (s - a) 2 ] - 38(6 - a)(s - 6) 2 + -(6 - af 



ds. 



Let us denote the integrands of the right-hand side of the above relation by 
<3i(s), Q2{s), Qs(s) and Q^(s), respectively. These integrands have the next 
zero points 



s = a, Si 



38a + 76 2a + 6 a + 26 



45 



S3 



-, S 4 



7a + 386 
45 ' 



sk = 6. 
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From the above two relations we get 

b 



! 



7g(o, 8) + 32q(^±l, S ) + Uq(^, S ) + 32«(i±^, s) + 7q(b, s) 



ds 



J Q 1 {s)ds~ J Qi(s)ds- J Q 2 (s)ds+ / () 2 l,-)<l* 

a SI 3a + b 



Q 3 {s)ds 



3 (s)ds - / Q i (s)ds + / Qi(s)ds 



2 

2036 



6075 
Thus, we have 



{b-af 



|i?(a,6)|<^(r- 7 )(6-a) 3 . 



From EOJl and (J3TJI we easily get lt2^|) . 
Remark 9 27ie usual Peano error bound is 

[7/(o) + 32/(2^) + 12/(^) + 32/(^) + 7/(6)] 

6 



-Jf(t)dt 



< 



1018 II fll 



273375 WJ Moo 



iiriioo(6-or- 



(31) 



for i/ie reasons given in Remark\^the estimation obtained in Theorem^ which 
is a Peano-like bound, is better than the above Peano bound. 
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